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Abstract 

We study thermodynamic characteristica of ultra-relativistic Au-|-Au collisions 
at RHIC energy ,/sJfff = 200 GeV within the framework of a microscopic trans¬ 
port model, namely the quark-gluon string model (QGSM). The temporal evolution 
of the local energy density, transverse and longitudinal pressure and equilibration 
times are considered. In contrast to complete equilibration which is even in central 
reactions hardly achieved, pre-equilibrium stages with energy densities well above 
the critical energy density predicted by lattice QGD are established at short time 
scales. Corresponding energy density profiles at proper time r = 1 fm/c compare 
well with hydrodynamical assumptions for initial energy density distributions. 
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1 Introduction 


Ultra-relativistic heavy ion collisions as performed within experiments at the 
Relativistic Heavy Ion Collider (RHIC) in Brookhaven offer a nniqne oppor- 
tnnity to stndy the nnclear phase diagram at high temperatnres and densities 
[1]. For snch extreme conditions the fnndamental theory of strong interactions 
(Qnantnm Chromo Dynamics - QCD), i.e. lattice QCD calcnlations [2,3,4] 
predict a transition to a plasma of deconhned qnarks and glnons (QGP). This 
state is believed to behave like a strongly conpled liqnid rather than a free 
parton gas [5]. 

It is a generally accepted opinion that the preconditions for the phase transi¬ 
tion to a QGP are readily fnlhlled at RHIG energies. In the meantime exist 
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several experimental evidences that such a phase transition has indeed been 
seen [6,7,8,9]. 

The strongest theoretical argument comes from lattice QCD which predicts 
the phase transition at a critical temperature of around 170 MeV which corre¬ 
sponds to a critical energy density of around 0.7... 1.3 GeV/fm^ (depending 
on the number of degenerate quark flavours) at zero net baryon density [10,11]- 
Estimates from pseudorapidity distributions of charged particles and elliptic 
flow results from the PHOBOS Collaboration indicate that in central collisions 
at yCjvTv = 200 GeV energy densities e > 3 GeV/fm^ are produced [7] when 
the system reaches approximate equilibrium. Such high energy densities, which 
are about twenty times the energy density inside nuclei and about six times 
the energy density inside nucleons, exceed clearly the values for the critical 
energy density predicted from lattice QCD. However, lattice calculations are 
performed for equilibrated matter and inhnite systems. This rises questions 
for the size of the space-time volume in which a heavy ion collision exceeds the 
critical energy density and also the degree of equilibration which is reached 
inside this volume. Therefore the evolution and absolute magnitudes of the 
energy density within the hrst few fm/c after the collision are of particular 
interest. 

Also the fact that the system enters the hydrodynamical regime at RHIC 
energies yCiT/v = 130 -i- 200 GeV is often interpreted as an evidence for a 
phase transition [12]. In contrast to SPS, at RHIC a simultanous description 
of elliptic flow pattern as well as of low px (pt ^ 2 GeV) single particle spectra 
around mid-rapidity can be achieved within the hydrodynamics approach. 
This success has led to the conclusion that the created matter behaves like 
an ideal fluid with almost negligible viscosity [13,14] and that the time scale 
to reach such an locally equilibrated state is extremely short, i.e. less than 
1 fm/c [15]. 

From an experimental point of view, it is almost impossible to extract informa¬ 
tion about initial pressure or energy densities and their equilibration times. 
Hydrodynamical calculations themselves provide only little insight in these 
questions since they are based on the assumption of local thermal equilibrium 
and short equilibration times. However, also hydrodynamics does not provide 
a perfect description of RHIG data. In particular, deviations of the hydrody¬ 
namical predictions for V 2 and from the measured elliptic flow were recently 
interpreted as a breakdown of the hydrodynamical picture and a hint towards 
incomplete equilibration [16]. In the hydrodynamics approach an ansatz for 
the early stage, i.e. at about 1 fm/c after the initial collision, is used as input 
for the evolution of the system. This requires to fix the initial energy den¬ 
sity and net-baryon number density prohles. Usually for both distributions a 
boost-invariant distribution over a wide range of ps, the so called space-time 
rapidity is chosen. 
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Transport models on the other hand, allow to study the dynamical evolution of 
thermodynamical properties of the colliding system from a microscopic point 
of view. String-cascade models of that type, e.g. UrQMD [17,18], HSD [19,20] 
or the quark-gluon string model (QGSM), do not contain an explicit phase 
transition to a deconhned partonic state but they can be used to study the 
preconditions for a possible phase transition and the results can be compared 
with the assumptions made by hydrodynamical approaches. The quark-gluon 
string model, applied in the present work and described in detail in Refs. 
[21,22,23,24,25] is a string model based on Gribov Regge Theory. The main 
assumption within these kind of (string-) models is that hadrons are pro¬ 
duced as a result of excitation and decay of open strings with different quarks 
or diquarks at their ends. The QGSM provides a safe basis for the present 
investigations since it has been demonstrated to describe the elliptic flow 
bulk properties at SPS and in particular at RHIG energies reasonably well 
[26,27,28]. 


2 Evolution of energy density 


In order to obtain an impression of the temporal evolution of the energy 
density e we consider as two typical examples a central (5 = 0 fm) and a 
semi-peripheral (6 = 8 fm) Au-|-Au reaction at top RHIG energy of = 

200 GeV. In the present approach the energy density is locally determined 
in cells of 1 fm^ volume through the scalar mass density of all formed and 
pre-formed hadrons inside these cells ^. All QGSM results are obtained from 
averages over 200 central reactions (6 = 0 fm), respectively 600 semi-peripheral 
reactions (6 = 8 fm). For simplicity, kinetic contributions to the energy density 
are neglected. Thus, the e values presented in the following analysis have to 
be interpreted as lower limits of the total amount achieved in the collisions. 
To clarify this a bit more, let us consider a central collision at the very first 
instant when the two Lorentz contracted gold nuclei with a longitudinal size 
of ~ 0.1 fm and a transverse area of ~ 150 fm^ overlap totally. This yields an 
energy density of ~ 25 GeV/fm^ in this approach. 

Fig. 1 displays the time evolution of the local energy density e in the x-z plane 
as a function of the system time t in the c.m. frame of the colliding nuclei. 
One should thereby keep in mind that the corresponding hypersurfaces where 
e(x, f = const.) is evaluated are essentially different form hypersurfaces of 
proper time e(x, r = const.) which are used to characterize the time evolution 
in hydrodynamics. The latter case will be discussed later on. 

One immediately recognises in the contour plot of Fig. 1, that most particles 

^ Exceptions to this specification are appropriately indicated. 
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b = Ofm I =2GeV/fm^ b = 8 fm 



z [fm] z [fm] 


Fig. 1. Evolution of the energy density in the x-z plane for two different impact 
parameters, 6 = 0 fm (left) and 6 = 8 fm (right). Each line corresponds to an energy 
density interval of 2 GeV /ivo?. 

do not stay in the central collision zone but rather form some kind of “shock 
waves”, i.e. two disks in the fragmentation regions occupied with large numbers 
of produced particles which propagate with high rapidities in front of the 
expanding system. As a result of that, the energy density in these zones stays 
at a high level for the whole time under consideration and decreases only 
slowly. Therefore most collisions happen within these zones while the center 
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of the volume between the two disks is rapidly thinning out of hadrons. After 
the hrst initial collisions energy densities of a few tens of GeV are achieved in 
the overlap zone of the central as well as semi-peripheral reaction. 

To study the initial phase and the overlap zone more closely, we consider now 
the average energy density inside a larger volume of 64 fm^. For the central 
reaction we choose the central reaction zone, i.e. the most central cell given 
by —2 fm < X, y, 2 ; < 2 fm and for the 6 = 8 fm case a volume shifted by 
4 fm in x-direction, i.e. 2 fm < x < 6 fm, —2 im. < y,z < 2 fm. Fig. 2 dis¬ 
plays the time evolution of e in the corresponding cells. One sees, that for the 
central collision a maximum energy density of around 15 GeV/fm^ is reached 
while the semi-peripheral collision reaches a maximal e of around 10 GeV/fm^. 
Again, we want to emphasize that all particles are taken into account, i.e. hrst 
of all no rapidity cut is applied. The energy density stays almost constant 



tvol 


Fig. 2. Evolution of the energy density for a central cell in the overlap zone of the 
collision given by —2 fm < x,y,z < 2 fm for 6 = 0 fm and 2 fm < x < 6 fm, 
—2fm<y,z<2fm for 6 = 8 fm. The curves with open symbols are the same but 
only for particles with a pseudorapidity of |? 7 | < 1. See discussion in the text. 

during the hrst 2 fm/c. It is even slightly increasing due to the ongoing par¬ 
ticle production through scattering of projectile and target nucleons and the 
corresponding continous energy deposit. As soon as the fastest particles start 
to leave the cell, what happens at around t ^ 2 fm/c, e is rapidly decreasing. 
This means that we obtain energy densities of 15 GeV/fm^ in a space-time 
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volume of at least ~ 130 fm^ in the central collision and correspondingly about 
10 GeV/fm^ in the same volume for the semi-peripheral reaction. However, 
at that time the system is far from local equilibrium, as will be seen later 
on. It is stated in Ref. [7] that the value of the energy density for such a 
system, although well dehned, may not be very interesting. The potentially 
more interesting quantity should be the energy density carried by particles 
which are closer to equilibrium conditions, wherefore the analysis in [7] was 
restricted to hadrons with pseudorapidities |? 7 | < 1. To make a comparison 
between the estimates of the PHOBOS Collaboration and the QGSM results 
more transparent, two additional curves for the temporal evolution of e with 
the same |? 7 | < 1 restriction are depicted in Fig. 2. Here, energy densities in 
the order of > 4 GeV fm“^ in the central collisions and > 2 GeV fm“^ in the 
semi-peripheral reactions are obtained after 2 fm/c time of evolution. These 
values are in nice agreement with the lower limit estimation of e > 3 GeV fm“^ 
produced in central Au-I-An collisions at RHIC after 1 — 2 fm/c, as it is re¬ 
ported by the PHOBOS Collaboration [7]. Note, that also the QGSM results 
are lower limits of the produced energy densities as mentioned before. The 
actual densities could easily be signihcantly larger. Nevertheless, the parton- 
string cascade model is not able to equilibrate completely the system at such 
short time scales, as we will discuss later on. 

As clearly seen from Fig. 1, the geometry of the colliding system is dominated 
by the leading particles which move practically on the light cone in beam [z) 
direction. In order to remove distortions from Lorentz contraction effects in a 
rapidly longitudinally expanding system one can introduce the so called space- 
time rapidity r]s. This is done by the introduction of cylindrical coordinates in 
t- and 2 ;-direction according to 

{t,x,y,z) = {TCOshr]s,x,y,Tsmhr]s) (1) 


with the longitudinal proper time r = — z"^. The space-time rapidity is 

then given by the following expression: 

1 t -\- z 

r]s = - In - 

' 2 t-z 

The evolution of the energy density can now be displayed as a fnnction of 
the boost-invariant variable r]s- Fig- 3 shows the energy density in the x — rjs 
plane at the same time steps as have been used in Fig. 1. The y-coordinate 
is thereby restricted to the interval 0 < y < 1 fm. Also as a function of rjs 
the system has an energy density distribntion with a nontrivial shape in the 
reaction plane. For the central collision, e is rather Gaussian distributed in 
x-direction whereas it shows a characteristic donble-peaked structure in rjg- 
direction, where the maxima follow the leading particles. This qualitative fea¬ 
ture of the energy density distribution is practically maintained in the overlap 
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Au+Au, s'^^=200 GeV, 0<y<1 fm, b=0 fm, t=1 fm/c 
E [GeV/fm^] 



Au+Au, s'^^=200 GeV, 0<y<1 fm, b=8 fm, t=1 fm/c 
e [GeV/fm^] 



Au+Au, s^^^=200 GeV, 0<y<1 fm, b=0 fm, t=4 fm/c 
E [GeV/fm^] 



Au+Au, s'^^=200 GeV, 0<y<1 fm, b=8 fm, t=4 fm/c 
E [GeV/fm^] 




Au+Au, s^^^=200 GeV, 0<y<1 fm, b=0 fm, t=8 fm/c 
E [GeV/fm^] 

2 r 



Au+Au, s^^^=200 GeV, 0<y<1 fm, b=8 fm, t=8 fm/c 
E [GeV/fm^] 

1.5 r 


Au+Au, s'^^=200 GeV, 0<y<1 fm, b=0 fm, t=12 fm/c 
E [GeV/fm^] 



1.25 



Au+Au, s'^^=200 GeV, 0<y<1 fm, b=8 fm, t=12 fm/c 
E [GeV/fm^] 

1 r 

0.75 


Fig. 3. Snapshots of the energy density e in the x — r]s plane at times 
t = 1, 4, 8, 12 fm/c. Note the different scales for e. 

zone of the semi-peripheral collision, albeit it is narrower distributed along the 
x-axis. Additionally, one clearly recognises the e contributions of the spectators 
(“shoulders” at large |x| and Ir^^l) which propagate with about the same veloc¬ 
ity as the high density bumps formed by the leading particles. Within the hrst 
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few fm/c of the expansion the central valnes are large, i.e. e ~ lO-i-15 GeV/fm^ 
and the distribntions are more compact. However, as already seen from Fig. 
2, e drops rapidly when the system expands. This decrease with increasing 
time is more prononnced in the central zone forming the “dip” than for the 
bnmps following the leading particles. Bnt also there the absolnte valnes drop 
signihcantly and the peaks get washed ont dne to ongoing secondary collisions 
between the prodnced hadrons. 

Complementary to the energy density in the reaction plane, the corresponding 
distribntions of e in the transversal (x-y) plane are depicted in Fig. 4. Here 
we consider only the central reaction with 6 = 0 fm. The contribntions of the 


Au+Au, s'^^=200 GeV, |t 15|<0.125, b=0 fm, t=1 fm/c Au+Au, s'®=200 GeV, |t 13|<0.125, b=0 fm, t=4 fm/c 

E [GeV/fm^] e [GeV/fm^] 





Au+Au, s’'^=200 GeV, |t 13|<0.125, b=0 fm, t=8 fm/c 
E [GeV/fm^] 

1 r 
0.8 - 


Au+Au, s’'^=200 GeV, |t 13|<0.125, b=0 fm, t=12 fm/c 
E [GeV/fm^] 

0.3 r 


Fig. 4. Snapshots of the energy density in the transversal plane for a space-time 
rapidity slice of \r]s\ < 0.125 obtained for a central (5 = 0 fm) Au-|-Au collision at 
times t = 1, 4, 8, and 12 fm/c. 

leading particles are snbtracted by restricting the space-time rapidity to the 
interval Irjsl < 0.125 in all the fonr time steps nnder consideration. The distri¬ 
bntions look now more isotropic, i.e. Ganssian shaped in x- and y-direction, 
and are peaked at the origin of the transversal plane. As already seen from the 
previons hgnres, the peak valne of e reaches np to 15 GeV/fm^ for t = 1 fm/c. 
Althongh the energy density is rapidly decreasing in the transverse directions 
as well as in time, there exist dehnitely regions with valnes of e within the hrst 
8 fm/c, that are higher than the expected critical energy density for a phase 
transition to a deconhned phase. 

Finally, the energy density e is evalnated for the central reaction at constant 










proper time r = 1 fm/c. The QGSM result is shown in Fig. 5. The depicted 
distributions can be directly compared with initial energy density distributions 
used in hydrodynamical calculations [29,30,31]. The e distribution in the reac- 


Au+Au, s'^^=200 GeV, |y|<1 fm, b=0 fm, t=1 fm/c 
e [GeV/fm^] 



Au+Au, s^^^=200 GeV, |Ti5|<0.125, b=0 fm, t=1 fm/c 
s [GeV/fm^] 



Fig. 5. Energy density of central Au+Au collisions in the reaction plane (|y| < 1 fm 
and \x\, \r]s\ > 0) ~ left panel, and in the transversal plane (|r/s| < 0.125 
and |x|, |y| > 0) - right panel, evaluated for constant longitudinal proper time 
r = 1 fm/c. 

tion plane as obtained in the QGSM is rather similar to the initial energy den¬ 
sities investigated by Hirano et al. for Au+Au collisions at ^/s]^ = 130 GeV 
in a hydrodynamical model [29,30]. The distribution is nearly flat in the re¬ 
gion |? 7 s| < 1.0 and connects to vacuum as a kind of Gaussian function sharply 
decreasing in the forward and backward space-time rapidity regions. The di- 



Fig. 6. Energy density in the reaction plane of central collisions obtained in the 
QGSM for T = 1 fm/c - left panel (the same as on the left of Eig. 5 but rescaled), 
and initial distribution used in hydrodynamical calculations by Hirano et al. - right 
panel. The figure on the right is taken from Ref. [29]. 

rect comparison of the microscopically determined density prohle on the one 
side and a hydrodynamical model assumption from Ref. [29] on the other side 
- illustrated in Fig. 6 - exemplihes this fact. Also in the transversal plane 
the shape and order of magnitude of e at r = 1 fm/c is supported by a hy¬ 
drodynamical calculation of Ref. [31], where the NeXuS event generator [32] 
has been used to determine the initial conditions for central Au+Au colli¬ 
sions at ^/SNN = 130 GeV. There it has been shown that fluctuations of the 
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initial conditions produced by a microscopic model and used as input for hy- 
drodynamical calculations are important and lead to discrepancies for some 
observables. Together with our hndings for the averaged distributions of en¬ 
ergy density the assumption of smooth and symmetric initial conditions seems 
to be a rather rough one. 


3 Local equilibration 


As one can expect from the energy density distributions of Fig. 1, global 
equilibrium cannot be reached; the distributions of particles are far from being 
isotropic for all time steps under consideration. But when one looks in more 
detail, some cells could reach local equilibrium. In order to measure the degree 
of kinetic equilibration for a given cell, the pressure components for the x-, 
y- and z-direction are calculated [33]. The diagonal elements of the pressure 
tensor in a cell according to the virial theorem are given by [34] 




{x,y,z} 


1 


E 




raf + 


l,X 


+ ply + P\ 


i,z 


( 3 ) 


where V is the volume of the cell, rrii is the mass and Pi{x,y,z} are the momentum 
components of hadron i . The sum runs over all hadrons in the correspond¬ 
ing cell. The pressure is locally evaluated, therefore a Lorentz transformation 
with the average velocity in the x-, y- and z-direction is applied before taking 
the sum. The requirement of isotropic velocity distributions, i.e. kinetic equi¬ 
librium, is closely related to the requirement of the pressure isotropy, which 
would imply = Py = Pz ■ Therefore, the equilibration ratio i?LE as dehned 
by 


.Rle 


1 Px P Py 

2 Pz 


( 4 ) 


is evaluated for each cell as a measure of local equilibration. Local equilibrium 
in the corresponding cell is then characterized by i?LE —^ 1 • 

First of all, the equilibration process is studied for the central cell of the col¬ 
lision, given by —2 im. < z < 2 im. for & = 0 fm and 2 fm < a; < 6 fm, 
—2 im. < y,z < 2 im. for b = 8 fm, by looking at the pressure for all three 
directions. Fig. 7 shows the pressure components as dehned in Eq. (3) for the 
central and the semi-peripheral collisions. In both reactions the longitudinal 
pressure is clearly dominating the very early evolution of the system while 
the transversal pressure components in x- and y-direction are nearly equal. 
At f ~ 2 fm/c the longitudinal pressure starts to decrease drastically and 
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Fig. 7. Pressure components in x-, y- and z-direction in the central cell of the overlap 
zone of the colliding gold nuclei with 6 = 0 fm and 6 = 8 fm, respectively. 

reaches a level close to the transversal pressure after 4 fm/c time of evo¬ 
lution. Thereafter all the three corresponding pressure components decrease 
with similar slopes. Nevertheless, the pressure isotropy is not yet achieved 
because deviations up to 30 — 50 % between the longitudinal and transversal 
components occur during the next 8 fm/c . Following the discussion in [33], 
one can call this a pre-equilibrium stage . It sets in very quickly for both, the 
central and the semi-peripheral collision, whereas the dynamics of the kinetic 
equilibration process during the pre-equilibrium stage is somewhat different 
for both reactions. 

This difference in the degree of equilibration is clearly seen, when one addition¬ 
ally looks at the equilibration ratio i?LE dehned in Eq. 4. Fig. 8 depicts i?LE 
in the central cell of the overlap zone for central and semi-peripheral And-An 
collisions. Again, the pre-equilibrium stage with R^e ~ 1 is reached for central 
collisions after 4 fm/c corresponding to a time when the high density areas 
as seen in Fig. 1 and Fig. 2 have completely passed through the central cell. 
Within these high density areas the longitudinal pressure is clearly dominat¬ 
ing and therefore no equilibrium can be reached. During the next few fm/c 
the transverse pressure components are slightly dominating the expansion; the 
equilibration ratio tends to be greater than one. But it decreases already after 
6 fm/c time of evolution and - as the cell is loosing more and more particles 
- approaches almost local equilibrium after 12 — 15 fm/c. The transverse pres- 
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Fig. 8. Degree of equilibration as defined by the equilibration ratio i?LE in the 
central cell of the overlap zone for Au+Au collisions at ^sjsfisr = 200 GeV with 
impact parameters 6 = 0 fm and 6 = 8 fm. The triangles indicate the equilibration 
ratio of particles with pseudorapidities |r/| < 1 only, produced within the first 3 fm/c 
of the central collisions (see detailed discussion below). 

sure components Px and Py are nearly equal to each other during this stage, 
what reflects the symmetry of the collision zone in the transverse plane. 

In case of the semi-peripheral collision with an impact parameter of 8 fm local 
equilibrium seems not to be reached in the cell corresponding to the overlap 
zone of the collision after 15 fm/c, but nevertheless i?LE gets at least very close 
to unity. After the high density areas passed through the cell, the transverse 
pressure is clearly dominating the next few fm/c of the pre-equilibrium stage 
- stronger than in the central collision, but in contrast to the central reaction 
it is itself not equilibrated. As can be seen in Fig. 7, the y-component of 
the pressure tends to be greater (up to 30%) than the x-component . This 
asymmetry arises due to the elliptical overlap of the colliding nuclei and reflects 
the smaller pressure gradients in the y-direction, which, in the hydrodynamical 
picture, is responsible for the build-up of elliptic flow [35]. 

The question, to what extent the produced matter reaches thermalization and 
how fast this may be achieved, is still open. In Ref. [16] it was shown that 
the linear rise of elliptic flow, in particular V 2 /e with e being the eccentricity 
of the spatial overlap of the colliding nuclei at RHIC, with the density gives 
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rise to the conclusion of a density dependence and hence to a small number of 
collisions per particle. Thermalization could therefore not be reached so fast. 
Our hndings of equilibration at around t ~ 15 fm/c together with the elliptic 
flow results and their time evolution obtained in the QGSM [28,36,37] support 
this picture and contradict quasi the usual hydrodynamical assumptions of fast 
thermal equilibration. 

However, as can be seen in Fig. 8 a certain fraction of the produced particles, 
i.e. hadrons with pseudorapidities |? 7 | < 1, is closer to local equilibrium even 
at times t ^ 1 fm/c. In between the two currents of matter of the target and 
projectile fragmentation regions, which make up ~ 2/3 of the energy density, 
this fraction provides a dense and roughly equilibrated local surrounding even 
at very early times of central heavy ion collisions. Thus we hud energy densities 
of about 4 — 5 GeV fm“^ in central Au+Au reactions at yGjvTv = 200 GeV 
for matter which can be considered to be at least in a pre-equilibrium stage 
already at 1 — 2 fm/c. The preconditions for a phase transition, i.e. energy 
densities e > Ccrit — 1 GeV fm~^, are given for a duration of ~ 8 fm/c in a 
space-time volume of about 500 fm"^. 

Finally, the degree of local equilibration and its time evolution is evaluated in 
the whole reaction {x — z) plane for the central and semi-peripheral Au-|-Au 
collisions. All the results depicted in Fig. 9 for the three time steps t = 
4, 8, and 12 fm/c show a characteristic saddle structure, which is more 
smoothly and faster developed in the central than in the semi-peripheral re¬ 
action. For the x-direction the equilibration ratio R^e always has a minimum 
in the overlap zone whereas it reaches here the maximum as a function of the 
longitudinal coordinate z. The high density areas in front of the expanding 
system, in other words at large values of z, never reach local equilibrium; only 
in the aftermath of the pass-through the inner cells tend towards kinetic equi¬ 
libration indicated by Rle ~ 1. In the center of the collision the transverse 
pressure is dominating at all times, therefore is > 1. But towards higher 
densities in the z-direction the longitudinal, i.e. z-component of the pressure 
is increasing and therefore supplying more equilibrated zones hlled with out¬ 
ward moving particles. Since the most collisions naturally happen within the 
high density areas, it is important to notice that the most secondary particle 
interactions take place in a non-equilibrated surrounding. 

Glosing the discussion, it should be mentioned that string-cascade models pro¬ 
vide most likely an upper limit for relaxation times towards local equilibrium. 
This is not only due to the fact that models of this type do not contain a 
phase transition to partonic degrees of freedom which are expected to equi¬ 
librate much faster. Standard string-cascade models violate detailed balance 
since particle production through string decays leads to processes 2 —N 
with two initial and N hnal state hadrons. However, corresponding backward 
processes N —2 are not included. This approximation is usually justihed by 
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Fig. 9. Snapshots of the equilibration ratio Rle in the overall x-z plane after 
t = 4, 8, and 12 fm/c time of evolution. The distributions on the left are for 
collisions with 6 = 0 fm, and on the right for collisions with 6 = 8 fm. 

the assumption that many-body scattering processes N —>■ X play a minor 
role in open, fast expanding systems. However, when high particle densities are 
achieved such processes can become relevant, e.g. for the production of rare 
hadron species [38]. Many-body processes, such as string or hadron fusion and 
recombination, lead naturally to shorter equilibration times, as has e.g. be 
demonstrated in [39] for the gluonic sector. In very dense states of matter it is 
expected that interactions among strings, before they fragment into hadrons, 
are able to accelerate the equilibration process in microscopic string models. 
Possible concepts in this context are the before mentioned string fusion, the 
formation of string clusters, and eventually percolation of strings, which have 
been considered in several string-cascade models [40,41,42,43,44]. Effectively, 
such processes describe reactions of 2n {n > 2) initial particles going to N 
hnal state hadrons. Nevertheless, detailed balance is still violated. But fu- 
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ture studies should be dedicated to the question how far these multi-string 
processes are able to reduce quantitatively the relaxation times towards local 
equilibrium. 


4 Summary and conclusions 


In this survey, we have analysed several aspects of the collision dynamics in 
central and semi-peripheral Au-|-Au reactions at the highest RHIC center of 
mass energy of = 200 GeV within the microscopic quark-gluon string 

model (QGSM), which is based on the colour exchange mechanism for string 
formation. 

The first major topic was dedicated to the time evolution of the local energy 
density e. Maximum energy densities of around 15 GeV/fm^ for the central 
collisions and around 10 GeV/fm^ in semi-peripheral reactions are reached in 
the overlap zone after the first initial hadron collisions. The energy density 
stays almost constant in the central cell during the hrst 2 fm/c time of evo¬ 
lution. The value of e in this cell is rapidly decreasing as soon as the fastest 
particles start to leave the cell. This happens at around t ^ 2 fm/c. The cor¬ 
responding e distributions in the x-z plane as a function of the system time t 
in the c.m. frame of the colliding nuclei reveal, that most particles do not stay 
in the central collision zone but form two disks in the fragmentation regions 
occupied with large numbers of produced particles which travel practically on 
the light cone in beam [z) direction. Also as a function of the boost-invariant 
space-time rapidity r]s, the system has an energy density distribution with 
a nontrivial shape in the reaction plane, i.e. a characteristic double-peaked 
structure in ? 7 s-direction which is maintained in the central and also in the 
semi-peripheral collisions. Gomplementary to the energy density in the reac¬ 
tion plane, the corresponding distributions of e in the transversal plane and 
additionally at constant proper time r = 1 fm/c were studied (only for the 
central reaction). The QGSM results of the latter distributions compare well 
with hydrodynamical assumptions for initial energy density prohles [29,30,31]. 
These are similar in shape as well as in the order of magnitude. Nevertheless, 
the absolute values of e drop rapidly when the system expands and the peaks 
get significantly washed out due to ongoing secondary collisions between the 
produced particles. But we want to emphasize that, although e is rapidly de¬ 
creasing in the transverse directions as well as in time, there dehnitely exist 
regions with local energy densities higher than the expected critical energy 
density for the QGD phase transition within the hrst 8 fm/c. 

However, at that time the system is not completely in local equilibrium, as 
it was elaborated in the second major topic. There the question whether and 
how fast the produced particles in ultra-relativistic Au-|-Au collisions equili- 
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brate has been studied. It was shown that in both reactions, i.e. central and 
semi-peripheral ones, the longitudinal pressure is clearly dominating within 
the first 2 fm/c. After the fast particles have left the central cell, all the three 
pressure components converge to each other. The next 8 fm/c the transverse 
pressure is slightly stronger, not allowing the system to reach full local equi¬ 
librium but entering some kind of pre-equilibrium stage. As the central cell of 
the overlap zone is loosing more and more particles, it approaches almost local 
equilibrium after 12 — 15 fm/c in the central collisions and 15 — 20 fm/c in 
case of the semi-peripheral reactions. Finally, the degree of local equilibration 
characterized by Rle and its time evolution in the whole reaction [x — z) 
plane shows in the central and semi-peripheral Au-|-Au collisions for all con¬ 
sidered time steps a characteristic saddle structure, which indicates that the 
high density areas in front of the expanding system are far away from local 
equilibrium. Only the inner cells tend towards kinetic equilibration. Since the 
most secondary collisions happen within these high density areas, the most 
particle interactions take place in a non-equilibrated surrounding. Global equi¬ 
librium is not reached at all - the global particle distributions are far from 
being isotropic for all time steps under consideration. 


As it was already mentioned, the question to what extent the produced matter 
reaches thermalization and how fast this may be achieved in the experiment 
is still open. Our hndings of almost full local equilibration at around t ~ 
15 fm/c together with the elliptic flow results and their time evolution obtained 
in the QGSM [28,36,37] support the idea that fast thermal equilibration is 
not necessary to produce large elliptic flow, what contradicts somehow the 
usual hydrodynamical assumptions. Nevertheless, a pre-equilibrium stage with 
energy densities well above the critical energy density is already achieved at 
early time scales. This means that the preconditions for a phase transition can 
be found within the string-cascade approach which, due to the lack of multi¬ 
string and/or multi-particle fusion processes, provides an upper limit for the 
equilibration times. 


In this context, it is also important to stress, that all calculations are per¬ 
formed by averaging over 200 and 600 events for central and semi-peripheral 
reactions respectively. For a single collision event the energy density as well 
as local equilibrium is subject to strong fluctuations. It is not a priori clear 
to what extent these fluctuations may contribute to the hnal observables, as 
has been pointed out in Ref. [31]. This question together with a more detailed 
understanding of the freeze-out scenario for the hadron production and the 
influence of multi-string/particle reactions on equilibration deserves further 
investigations. 
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